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Abstract 

Generalized log-sine functions Lsj(0) appear in higher order e-expansion of differ- 
ent Feynman diagrams. We present an algorithm for the numerical evaluation of these 
functions for real arguments. This algorithm is implemented as a C++ library with 
arbitrary-precision arithmetics for integer < k < 9 and j > 2. Some new relations 
and representations of the generalized log-sine functions are given. 
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PROGRAM SUMMARY 

Title of program: Isjk 

Version: 1.0.0 Release: 1.0.0 Catalogue number: 

Program obtained from: http://thsunl.jinr.ru/~varg/dist/ 



E-mail: varg@thsunl.jinr.ru 

Licensing terms: GNU General Public Licence 

Computers: all 

Operating systems: POSIX 

Programming language: C++ 

Memory required to execute: Depending on the complexity of the problem, at least 32Mb 
RAM recommended. 

Other programs called: The CLN library for arbitrary-precision arithmetics is required at 
version 1.1.5 or greater. 

External files needed: none 

Keywords: Generalized log-sine functions, Feynman integrals, Polylogarithms. 

Nature of the physical problem: Numerical evaluation of the generalized log-sine functions 
for real argument in the region < 6 < tt. These functions appear in Feynman integrals. 

Method of solution: Series representation for the real argument in the region < 6 < tt. 

Restriction on the complexity of the problem: Limited up to ~Ls*f\6), and j is an arbitrary 
integer number. Thus, all function up to the weight 12 in the region < 6 < tt can 
be evaluated. The algorithm can be extended up to higher values of k (k > 9) without 
modification. 

Typical running time: Depending on the complexity of problem. See Table ^ in section 13.41 
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LONG WRITE-UP 



1 Introduction 



Quite recently it has been realized that generalized log-sine functions [1] play an impor- 
tant role in analytical calculations of multiloop Feynman diagrams. Within dimensional 
regularization [2] with an arbitrary space-time dimension d = 4 — 2e, the log-sine functions 
Lsj(#) appear in the construction of the all order e-expansion of the one-loop propagator 
type diagrams with arbitrary masses [3], three-point integrals with massless internal lines 
and arbitrary (off-shell) external momenta [4], two-loop vacuum diagrams with arbitrary 
masses [4,5] and two-loop single scale diagrams of propagator type [6]. The appearance of 
the generalized log-sine functions Lsf } (0) (k = 1,2) was first detected in the finite part of 
the three-loop bubble master- integrals [7,8]. More examples of the one-, two- and three- 
loop Feynman diagrams whose e-expansion contains the generalized log-sine functions can 
be found in [8-13]. The generalized log-sine functions are also related to the derivatives of 
the generalized hypergeometric functions with respect to their parameters [13,14] and to the 
multiple Euler-Zagier sums [7, 8] . A collection of the known analytical properties of these 
functions can be found in Appendix A of [8] (see also [15]). 

Here we present C++ code for arbitrary precision numeric evaluation of the generalized 
function Lsj (9) of the real argument < 6 < tt and k < 9. For calculation of the Ls 3 -(0) 

and hsj (9) functions there exists a FORTRAN program (see the description in [16]). 

The plan of the paper is the following. Section [2] contains the definition and some 
properties of the generalized log-sine functions. Section El is the description of the Isjk 
library. In Appendix El we present relations between the functions Lsj 1 (6) and infinite 
series containing the trigonometric functions and harmonic sums. 

2 Generalized log-sine functions 

2.1 Definition and reflection symmetries 

The generalized log-sine functions are defined as [1] 



o 

where k, j are integer numbers, k > and j > k + 1 and 9 is an arbitrary real number. 
According to the definition (|2.1|) . the following integral and differential relations hold: 




(2.1) 




(2.2) 



o 



^s^\a9) = (a9) k ^Lsf\a9) , 



(2.3) 
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where a is a real number. In particular, 

Ls?" 1 ^) = ■ 
J 

It is easy to obtain a relation between functions of the opposite arguments 

Lsf \-9) = (-l) k+1 Lsf\6) , 6>0. (2.4) 
Another relation was deduced in Ref. [1] via the integral representation: 

l4\2m«-0) = LsM(2m7r) + (-l)- 1 LsM(^)-E(-l) r ^(2m7r) p ( r ) Ls£j> (9) , (2.5) 



P =i 



where m is an integer number. Using the symmetry property (|2.4|) . this relation can be 
written as 

L S M(2m7r+0) = LsM(2m7r) + LsW(^) + E(-l) r ^(2m7r) p ( r )l^j\e) . (2.6) 



P =i 



V) 



The generalized log-sine functions appear in the decomposition of polylogarithms Lij (l — e ie ^j 
into real and imaginary parts (see Eqs. (A. 7) in [8]). It is easy to show that the generalized 
Nielsen polylogarithms [17] S a ,b{z) of a unit-circle complex argument z = e ld also reduce to 
a combination of the generalized log-sine integrals (|2.1j) 1 

M e *) = /-~1VW f** V-ti*' 1 [ln|2sinf|-|i(7r-0)]Vx:^r^(l) 



fc=0 



+ E^f <Wl) • (2-7) 
fc=o K - 

2.2 Identities between the generalized log-sine and Clausen func- 
tions 

It was shown in [8] that the function Ls£_ 2 (^) is always expressible in terms of the Clausen 
functions CL, (6). This follows from the fact that Ls2(#) = CI2 (9), the differential identity 
(|2.3|) . which in this case has the following form: 

^LsK(0) = 0^Cl 2 (0) , (2.8) 



lr The corresponding equation in Ref. [10] contains an misprint in the last term. We are grateful to 
A. Kotikov for correspondence. 
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and the integration rules for the Clausen function [1], 



Cl 2n (9) = J &<p Cl w (0) , Cl 2n+1 (9) = C 2 „+i - J d0 Cl 2n (0) 
o o 



The general solution of Eq. ()2.8j) could be represented in the following form 



-iy(2j + l)![Cl 2 , +3 W-C2 J+ 3] 
fc=0 l / J + i_/c J ! 



The corresponding expressions up to Lsg 4 ' ) (#) have been presented in [8] . 



(2.9) 



(2.10) 



2.3 Identities between the generalized log-sine functions and the 
generalized Nielsen polylogarithms 

In Ref. [8] the following relations between the log-sine functions and generalized Nielsen 
polylogarithms were deduced: 



ia [Lsj(7r) — Ls^ 



$3 



1 



■IV 



- 1)! E -^rr E (-2)- fc [Skj-k-vW - {-iys k ^ p {i/ z )\ , (2.11) 

p=0 / P- k=l 



where the relation between z and 9 is defined via 

z = e iae , ln(-z - iaO) = \n(z) - ian, a = ±1, 



(2.12) 



and < 9 < it. For higher values of j the number of generalized polylogarithms involved in 
Eq. (I2~TT1) can be reduced (see Eqs. (2.21), (2.22) in [8]). The relation between Ls£l(0) and 
the generalized Nielsen polylogarithms is also known (see Eq. (A. 20) in [8]), 



LsgiW-Ls^iW = ^ 



(i) 



Lsj(0) - 



\n j+1 (-z) 



2^'0' + l) 



1 



i~ 2 ] n P(- Z ) J '~ 1_p (-l) k 



p=0 

J'-l 



+2(-l)'-(j-l)!£ 



2Pp! 
f-lV 



fc=i 



-JfeS. 



fe+i,j-M 



-1). 



(2.13) 



fc=i 
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Using the algorithm described in [8] we found the following relation for hs j^(9) and gener- 
alized Nielsen polylogarithms: 



1(7 



Ls ? 2 (tt) - Ls fl 2 {9) 



2ia9 



1 



\n j+2 (-z) 



i - f-iV 



Lsj(vr) -Lsj(0) 



2i- 1 j(j' + l)(j + 2) 

-(-1)^-1)! ^ £ ±—L^ " [Xk~-2 :j -, 

P =o z ^- fe=l 



fc+2,j-fc 



(-1)' 



+4(-iy (j - 1)! ln(-z) £ K ~^kS k+1 



k=l 



't+l,j—k\ 



■I) . 



-p(iA)] 



(2.14) 



The procedure, described in [8], can be repeated several times, so that the generalized log-sine 
functions with arbitrary indices can be directly related to Nielsen polylogarithms. 



2.4 Series representation 

In the region < 6 < tt, the series representation of the generalized log-sine functions 
directly follows from the definition (j2.1j) by substituting y — sin | 



sin(0/2) 



Lsf (0) = -2 fc+1 / -=L= (arcsin^ln^- 1 ^). 

{ v 1 - r 

Using integration by parts, this integral can be rewritten in the following form: 



(2.15) 



sin(6»/2) 



V / y 



Expanding (arcsiny) fc in y (see below) and using 



x a \n b xdx=(-l) b b\x a+1 Y / ' h) '' 1 



b—p 



p f (6-p)!(a + l)^ ' 

we obtain the multiple series representation. The generating function for the coefficients of 
the Taylor expansion of (arcsin?/) fc and their explicit form for k = 1,2,3,4 can be found 
in [18]. We have calculated the coefficients of the Taylor expansion for k = 5, ■ • • , 12. For 
even powers, the Taylor expansion of (arcsin|/) 2fc can be represented in the following form: 



(arcsin y) 



2k 



g [(m ^^, )!]2 4 m -y m i?(///7,) 



(2*)! ^ (2m)! 
where the coefficients B(m, k) are expressible in terms of harmonic sums 

n ^ 

S M = E y« ■ 



k=l 



(2.17) 



(2.18) 



For the lowest values of k we have 

5(771,1) 
5(777,4) 

^(m, 5) 
5(m, 6) 



1, 5(m,2)=<7 1 , 5(m,3) = - 

1 r 3 1 

— \a\ - 3a t a 2 + 2a 3 J , 

1 



(T l - CT 2 



4! 
1 

5! 



+ 8(Ti(T3 + 3(T 2 — 604 — §(J x Oi 



o\ - 30o-i(T 4 - 2O020-3 + 20a? 0-3 - lOcrf a 2 + ISa^j + 24a 5 



and 



m-1 j j 
ff a = E 7777^ = — ^2a(r77 - 1) 



, - (2j) 2a 4« 

For odd powers, the Taylor expansion of (arcsin?/) 2fc+1 has the following structure: 



2fe+l 



(arcsiny) 
(2fc + l)! 



e r2ni 



m=0 



y 



2m+l 



m / 4 m (2m + 1) 



C(m,k + 1) 



where the coefficients C (m, k) also expressible in terms of harmonic sums, 



C(m,l) 
C(m,4) 
C(m, 5) 
C(m, 6) 



3! 
1 

4! 
1 

5! 



C(m, 2) =h 1: C(m,3) = - /7 2 - /7 2 



- 3/ll/la + 2/7; 



^ + 8/I1/I3 + 3/7 2 - 6/74 - 6tl% 



/if - 3O/71/74 - 20/72/73 + 20/7 2 /7 3 - 10/7? /i 2 + 15/ia/i^ + 24/i £ 



and 



m 1 1 
§ (2j - l) 2 " = 52a(2m ~ 1} ~ 4^ 2a( ™ " 1} ' 



Starting with Eq. (|2.16j) and using expansions (|2.17|) and (|2.21j) . we get 



r-2k-i r_ j n 

Lsf fc) (^) = (-l) J '(2sinf)2 2fc (2fc)!(j-2A;-l)! £ ± 

p=0 

^2m\ (2sin|) 2m 



2 sin 



pi 



00 

*E 

m=0 



-C(m,k + 1) 



Km J 16 m (2777 + iy- 2k -p 

fou^ J- 2fc -2 [- In f 2 sin I 

Lsf k+1) (6) = (-iy- 1 2 2k+2 (2k+l)\(j-2k-2)\ E 

1 (2sin|) 2m 



00 

xE 

m=l 



( 2 r) 



\j-2k-p- 



where k = 0,1,2,3,4,5 and the coefficients B(m,k) and C(m,k) are given in 1)2.19)1 and 
()2.22j) . respectively. 

The special interest in the calculation of the multiloop Feynman diagrams [7,8, 10, 19] is 
related to the value of the generalized log-sine functions at the "sixth root of unity", 6 = | . 

In this case In ^2 sin |) = and only the p = term in the finite sum in Eqs. ()2. 24)1 . ()2. 25)1 

survives, 

Lsf >(|) = (-l)^(2fc)!( J -2/fc-l)! £ ( 2m ) ,f^ k +t ] 2k , (2.26) 
L B f «'(|) = (-iy-'2»"(a+l)I(j-2*-2)l £ * ^|±i> . (2.27) 
Similar series for Ls.,(#) and Ls^\8) have been presented earlier in [16]. 



2.5 Special values of the argument 

The values of the generalized log-sine functions of the argument 9 = n and 2n can be 
related with the values of some well-known special functions. The relations of this type were 
discussed in Lewin's book [1] . For completeness we present some of them below. The values 
of Lsj(7r) can be expressed in terms of the Rieman (^-function, for any j (see Eqs. (7.108), 
(7.110) in [1]), 



Ls 



n+1 



7T 



— 7T 



dx 



r(i+x) 



r 2 (i-H 



x=0 



"Mdi exp 



.m=2 



[2.2%) 



J x=0 



so that 



Ls 2 (7r) = 0, Ls 3 (7r) = -\kC,2, Ls 4 (vt) = §7rC 3 , Ls 5 (tt) = -f ttC 4 , Ls 6 (tt) = f ttCs + f 71-C2C3, 



etc. Putting in Eq. (|2.5|) r = and 6 = 2ix{m — 1), solving the difference equation and taking 
into account the relation between values of function Lsj(#) at 6 = tt and 27r, we finally get 



From relation (|2.5|) for 



Ls n (2mvr) = 2mLs n (vr) . (2.29) 
1 and # = 27rm we have a new relation 

LsW(27rm) = 2m 2 7rLs„_ 1 (7r) . (2.30) 
Putting r = 3 and 9 = 2nm in Eq. (J23) and using Eqs. (jOH (l2~37IJ) we get 

Ls^(27rm) = 37rmLsi 2 2 1 (27rm) - 4m 4 7r 3 Ls n _ 3 (7r) , (2.31) 

where the value of Ls£f (2irm) should be calculated independently. The integral Ls^ (2tt) 
can be expressed in terms of a combination of (^-functions (see Eq. (7.136) in [1]): 



Ls 



(n) 



i(2tt) 



-27r(-ir ^ — 



d 

dx 



d 

Av, 



r(i+j/) 



r(l + |+a;)r(l + f-x) 



(2.32) 



2=^=0 
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It can be written in the following form 



Ls 



(n) 

n+m+1 



x exp 



oo yP 



^ + Ei[H) p (i-2i(p 

p=2 



g=2 



(2.33) 



- x=j/=0 



so that for n = 2 we get 



Lb£U(2tt) = -7r 2 Ls m+1 (7r)-7r(m + l)!(-2) 2 - m C 2+ r 



m—l 



y—2) P ~ 

+4m! V V- 



+m—p i 



m>2 



(2.34) 



In particular, we have 



i4 2) (2tt) = --7T 



13 



l4 2) (2tt) = 5tt 3 C 3 + 12ttC 5 , Ls^ (2 



(2), 



7T 



29 , 

7T 

105 



•24ttC 3 2 . (2.35) 



Substituting expression ()2.34|) in Eq. ()2.31|) we get (m = 1) 



Lsf(27r) = 47T 3 Ls n _ 3 (vr)-37r 2 (n-3)!(-2) 6 -"Cn-2 

+37r(n-4)!^^^ (n-p-3)Ls p+1 (7r) C„_ 2 - P , n>5, (2.36) 

P =i 



p! 



and 



Lsf (2vr) = 127T 2 C 3 , Ls^(2tt 



(3), 



vr 6 , LsS. 3) (2tt) = 9tt 4 C 3 + 36tt 2 C 5 • 



15 



(2.37) 



The values of the generalized log-sine functions of the argument tt are related to the 
values of the T-functions and generalized Nielsen polylogarithms of the argument z — — 1 
(see section 7.9.9 in [1]). For lower values of k (k = 0,1,2) the proper relations can be 
deduced from Eqs. fl£EQ ~ <EH : 



LSj(7r) 



2 j J 



~ 7F 



i-(-iy 



2 fc 



p=0 



5 fcii _ fe _ p (l) [1 - (-If] 



(2.38) 



fc=i 



2 The term -[l+(-l)"](n - 1)!C„ is absent in Eq. (7.139) of [1]. 

3 In Eq. (7.141) of [1], the term 27r 2 Ls r+ i(7r) should read §7r 2 Ls r+ i(7r). 

4 The correct result for LS5 (27r) was first given in [8]. 
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ii+i 



2^'(J + 1) 



7T 



fc=l z 



J '~ 2 fiTT^ (-1^ 

+(- w -i)!eS e [i + (-i)i 

p =o z k=l z 



(2.39) 



ii+3 



2^0' + 1)0' + 2) 



7T 



J+2 



3-1 



1_(_1)3 _4 7r( _ l) 3 (i _ l) ,£ 



fc=l 



2 fc 



k+l,j-k{ 



-1) 



3-2 

•iy(i-i)'E-r 

p=0 



7T 



) 3-1-P 

E • 

fc=i 



-l) fc £;(fc+l) 
2 1 



In particular, we have 

Lsl 1} (7r) = 
Ls^(Tr) = 



-25 2 , 2 (-l)-fC 



4 ■ 



65; 



2,31 



1) + |C 2 C3-#C 5 , Lsf (tt) = 27rLs^(7r)+27rC 4 



(2.40) 



(2.41) 
(2.42) 



where the values of generalized Nielsen polylogarithms of argument z = ±1 can be extracted 
from [17,20]. For completeness we present the proper results below 

SaM = <Mi) , S2, 2 (i) = id, 5 li3 (l)=C4, S 2>3 (l)=2{ 5 -{ 2 { 3 , 
S lt3 {-1) = ^ln 4 2-iC 2 ln 2 2 + |C3ln2-C4 + Li 4 (|) , 



(2.43) 



S2,s(-1) = -^ln 5 2 + |C 2 ln 3 2-|C3ln 2 2 + |C2C3 + lC5-2Li 4 (i)ln2-2Li 5 (i 
S 2)2 (-l) = 2S- 1)3 (-1) + |C4, 5-3,2(-l) = K3C2 - §C 5 • 



3 The Isjk library 
3.1 Domain 

Our algorithm for the numerical evaluation of the generalized log-sine functions in the region 
< 9 < 7r is based on the multiple series representation ()2.24|) and (|2.25j) . The explicit value 
of the generalized log-sine function at 9 = 71 and 9 = 2ir can be extracted by the approach 
described in Lewin's book [1] (see also section l2~o"j) . For an argument belonging to the region 
7r < 9 < 2ir the following representation is valid 

Lsf (9) = (-l) fc+1 E("2^) P Q {ls^vt " 9)-Lsfr;\2n)} . (3.1) 

For negative values of argument, 9 < 0, the reflection symmetry (|2.4|) is applied. Eqs. (|2.5|) 
and ()2.6|) allow us to express the value of the generalized log-sine function of the argument 
\9\ > 2tt in terms of the function Ls^ (0) at < 9 < n, ti and 2n. 

In the present version of Isjk we have considered only the region < 9 < ir. 
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3.2 Description 

For arbitrary-precision arithmetics, Isjk uses the CLN library 5 . CLN has a rich set of number 
classes: 

• Integers (with unlimited precision) 

• Rational numbers (with unlimited precision) 

• Floating-point numbers: 

— Short float 

— Single float 

— Double float 

— Long float numbers (with unlimited precision) 

• Complex numbers 

and implements many elementary functions on these numbers. 
Isjk provides two functions 6 

const cl_R Ls (const unsigned & j, const cl_R &x) ; 

const cl_R Ls(const unsigned & j, const unsigned & k, const cl_R &x) ; 

These functions represent Lsj(x) and LSj(x), respectively. The argument x is supposed to 
satisfy < x < tt, parameter k should be < k < 9 7 . Isjk can be used as any regular 
library, e.g. by writing programs and linking executables to it. 

3.3 Installation instructions 

The installation procedure is reduced to the simple three steps 8 of 
./configure && make && make install 

The "configure" script can be given a number of options to enable and disable various 
features. A few of the more important ones are documented in the "INSTALL" file from the 
Isjk distribution. 

After installation Isjk can be used as any regular library, e.g. by writing programs and 
linking executables with it. Isjk includes shell script "lsjk-config" that can be used for setting 
compiler and linker command-line options required to compile and link a program with the 
Isjk library. Thus, a simple program can be compiled as 

C++ -o simple 'lsjk-config — cppf lags' simple. cpp 'lsjk-config — libs' 
See "README" file from the Isjk distribution for more details. 

5 CLN can be downloaded fr om |http : / / www. ginac . de / C LN Alternatively, you may wish to check your 
favored operating system distribution for a precompiled package. 
6 cl_R is CLN type for real numbers, see CLN manual for details. 

'In the present version we have implemented only first five coefficients B(m, k) and C(m, k), k = 1,2, 3, 4, 5 
from Eqs. (|2.19|) and l|2.22[) . respectively. 
8 last step may require superuser privileges 
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J = 4 




3 = 


= 5 






3 = 


= 6 




d\k 





1 


2 





1 


2 


3 





1 


2 


3 


128 


0.053 


0.034 


0.055 


0.055 


0.037 


0.058 


0.048 


0.059 


0.04 


0.063 


0.052 


256 


0.166 


0.106 


0.187 


0.172 


0.111 


0.195 


0.16 


0.177 


0.118 


0.203 


0.166 


512 


0.685 


0.448 


0.872 


0.699 


0.47 


0.876 


0.742 


0.727 


0.491 


0.9 


0.741 


1024 


3.212 


2.219 


5.084 


3.294 


2.292 


4.459 


4.013 


3.327 


2.342 


4.496 


3.786 



Table 1: Dependence of the average calculation time (in second), as reported by getrusage(2), 
of Lsj fc) (27r/3) on the precision (d) of evaluation. It was measured on a Duron/800MHz with 
128Mb RAM. 



3.4 Benchmarks and cross-checks 

Correctness of implementation was cross-checked by comparison with results of numerical 
integration 9 of Eq. ()2.1j) for 4<j<ll,0<fc<9, with randomly chosen argument 
< x < 7r. We also checked that the results of the evaluation of the functions Ls j(6) and 
hSj(9) coincide with the corresponding results of the FORTRAN program described in [16]. 

To give a reader a rough idea about the evaluation time, Table[T]may be helpful. It should 
be noted that the calculation time for values of the argument near 7r increases substantially. 

3.5 Example of a program 

The following program evaluates Ls5^(27r/3) up to 1024 digits. 

#include <iostream> 
#include <cln/cln.h> 
#include <lsjk/lsjk.hpp> 
using namespace std; 
using namespace cln; 

int main(int argc, char** argv) 
{ 

const unsigned j=5; 
const unsigned k=2; 

f loat_f ormat_t prec=f loat_f ormat (1024) ; // precision 

cl_R ThePi = pi(prec); 

cl_R res = Ls(j, k, 2*ThePi/3) ; 

cout « res << endl; 

return 0; 

> 

The output of this program is: 



3 For numerical integration we used GNU scientific library, see http://www.gnu.org/software/gsl 
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-0.5181087868296801173472656387316967550218796682431532140673894724824\ 
6493059206791506817591796234263409228316887407062572713789701522832828\ 
8301238053344434601555482416349687142642605456956152340876808788330125\ 
2744524532005650653916633546607642565939433325023687049969640726184300\ 
7710801944912063838971724384311449565205834807350617442006639919209369\ 
6654189591396805453280242324416887003772837420727727140290932114280662\ 
5550331483934346157017999680014851656153847980079446225103428769237020\ 
8915162715787613074342208607956842722678511943161304156600252481686319\ 
3071905942805038395349632095408925090936865076348021184023670239583644\ 
8060572488328604823250125773056264305964195471500442032760480010663468\ 
6808425894080311707295955974893312124168605505481096916643118747707399\ 
7726995685152764342004787909957595709564369273416476440874928297302997\ 
3226232121625055566818301147295999943567467361944097333638320437023447\ 
2148571069312485138137784268429835976279026626925221417350916538992701\ 
34033573580034488504210827748452303627506494601119041L0 



4 Conclusion 

In this paper we have described the C++ library Isjk for arbitrary-precision numerical eval- 
uation of the generalized log-sine functions for the real argument < 9 < ir. The evaluation 
is based on the series representations 1)2.24)1 and 1)2.25)1 . The present version allows one to 
evaluate generalized log-sine functions up to Lsj(9), thus, all functions up to weight 11 can 
be evaluated. The algorithm can be extended to higher weights without modification. In 
particular, implementation of the next coefficients B(m, 6) and C(m, 6) (see Eqs. 1)2.19)1 and 
()2.22j0 allows one to evaluate all functions up to weight 14. 

Explicit formulae (|2.14)) relating Ls^ (9) and the generalized Nielsen polylogarithms have 

been obtained. Using the relations of this type for the functions Ls 3 (0) and Ls^(0), the 
relations (JA.2)) - (JA.4)) between series included trigonometric functions and harmonic sums 
and the functions Lsf' 1 ' 2 \9) have been deduced. As a particular case, some new identities 
(|A. 12)1 - (|A. 17)1 for a special type of harmonic sums have been obtained. Another application 
of these expressions is the relation between the values of functions Ls^- (9) at 9 = 7i and the 
Nielsen polylogarithms at z — ±1 (J2.39)) . ()2.40j) . The correct expressions for the functions 
Lsf' 3) (#) at 9 = 2tt in terms of the T-functions and their derivatives are obtained (j2.34j) . 

flUU). 
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A Identities between Ls»- ' (0) and infinite series 

Identities between the generalized log-sine functions and infinite series involving trigonomet- 
ric functions and harmonic sums follow from the power series expansion of the generalized 
Nielsen polylogarithms 



£>n,p( z ) — E 



k=l 



kW 



(A.l) 



where 



denotes Stirling numbers of first kind defined as [21] 



ln p (l - z) 
Using Eqs. (I2~TT|) and (IA~T| we get 



'WE 

k=p 



Z 



L Sj (9)-L Sj (7r) = —(-iy- 1 i 



2*j 

-(-1) J (J-1)!E 

p=0 



7T 



1 - (-1^ 



-10- 



71 



Z fc=l n=l 



j-k-p 



n\n" 



x {cos(n#) [l-(-l) p ]+iasin(^) [l + (-l) p ]}. 



(A.2) 



Similar expansion can be deduced for Lsj (0) and Ls? (0) with the help of Eqs. 1)2.13)1 and 

dm, 



1^(0) - Ls^tt) = 

-i-iyu-iy-T,- 

p=0 



,(1) 



Lsj (9) — Lsj (ir 



#j(j + l) 

p j—l—p ( 1 \k oo 



(tT-0) 



lV 

2*p! fe tl 2* ^ 



n 

j-k-p 



n\n 



k+l 



x {cos(n0) [l + (-l) p ]+iasin(n0) [l-(-l) p ]} 

+2(-iy'(j - l)!E^^ fc+1)i _ fc (-l) , 
fe=i z 



(A.3) 



Ls?i(0) - Ls^ 2 (tt) = 20 
p'+3 



,(2) 



Ls^-Ls^tt 



(i) 



-0 2 



Lsj(0) — Lsj(7r) 



71 



V'+2 



2^0' + 1)0' + 2) 
•( !)'(./ D!E 

p=0 



2Pp\ fe tl 2^ 

x {cos(n0) [l-(-l) p ] + io-sin(n0) [l + (-l) p ]} 



p j—l—p /_1 \fc oo 

- E 4r-^ + i)E 

fe=l z n=l 



n 

J-k-p 



k+2 



-1)' 



+4(tt - 0)(-iy (j - 1)! e 



fc=i 



(A.4) 
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For lower values of j these expressions have the following form: 



Ls 3 (0) = -2 £ ^i^Si - & (0 2 ~ 3*6 + 3tt 2 ) , 

1 



n=l 
oo 



Ls 4 (#) = 3 E 



smmfc 



n=l 



n- 



Sf-S 2 +Hir- 



\2 1 



2 n 2 



-|i7T-( 



v - COS(n„ 



n=l 



2^3 , 



(A.5) 
(A.6) 



Ls 5 (#) = 3(tt - 9) £ 



cosine 



n=l 



n- 



sin(ra6 l ) 



n=l 



rr 



-4E 



sm (nt 



n=l 



n- 



S 3 1 -3S 1 S 2 + 2S 3 \-±n9(n - 9)(n 2 -n9 + 9 2 )-±9 5 -±7i 5 
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Ls 4 » = -2]T 

n=l . 



cos(n0J ^sm^nfc 



rv 



rr 



51 - ^ 2 (6vr 2 - 8tt9 + 39 2 ) + ^tt 4 , 



(A.7) 
(A.8) 



ls«w = |E cosM)/6 8 



n=l 



— 5*1+4 S 2 — S 2 



n 4 n 



n=l 



2 ^2 



n 



S 2 -S 2 +IC2C3-IC5, (A.9) 



# 3 



n=l 



rr 



rr 



^2 CUS^ Si _ , (^2_ 15 ^ + 107r 2) (A1Q) 



where we have introduced a short notation S a = S a (n—1), and we used the general expression 
for the Stirling number of the first kind in terms of the harmonic sums [21] 



n 
.0. 


— Son , 


~n~ 


= J(n 


.3. 


~n~ 


= a(" 


.5. 


~n~ 




.6. 



(n- 1)! 



(n-1)! Sl-S 2 



--(n-l)\S 1 , 

f (n-l)! fe 3 - 3Si5 2 + 2S 2 



Using the results of paper [8] (see Eqs. (A.9), (A. 10), (A. 14) in [8]) we get from l)A.5| ) -(|A.9 |) 
the following relations: 



E 

n=l 
00 

E 

n=l 

00 

E 

71=1 



sin(| 




n 2 




cos (| 


-n) 


n 3 




sin(f 


n) 


n 2 



Si = ^(2 

C _ 23 / 
Dl ~~ 216 C ' 4 ' 



S 2 - S 2 



2C1 4 (f ) - |C 2 C1 2 (f ) - IttCs 



(A.12) 
(A.13) 
(A.14) 
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3E 



cos 



71=1 

g cos (f n) 
)i=i 



Sl-S 2 -2 



Si 



n 



Jn(fn) - N 

+ 27r E T-^-^i = 27rC1 4 (f) -IC2C3-IC5 , (A.15) 



XJ sm 

n=l 



rr 

M cos 



Si = -1C4 + 4 In 4 2 - A G ]n 2 2 + |Cs In 2 + §Li 4 (|) 



(A.16) 



3E 

71=1 



n" 



Sf-S 2 -2 



Si 



n 



E 

71=1 



sin (|n) 



Si - |7rCl 4 (f ) - y^Cs - §f C2C3 



-i In 5 2 + A C2 i n 3 2 _ i| C 3 ln 2 2 _ ^Li 4 (i) In 2 - f Li 5 (± 



(A.17) 
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